Fluctuation analysis of mechanochemical coupling depending on the type of 

bio-molecular motor 
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Mechanochemical coupling was studied for two different types of myosin motors in cells: myosin 
V, which carries cargo over long distances by as a single molecule; and myosin II, which gener- 
ates a contracting force in cooperation with other myosin II molecules. Both mean and variance 
of myosin V velocity at various [ATP] obeyed Michaelis-Menten mechanics, consistent with tight 
mechanochemical coupling. Myosin II, working in an ensemble, however, was explained by a loose 
coupling mechanism, generating variable step sizes depending on the ATP concentration and realiz- 
ing a much larger step (200 nm) per ATP hydrolysis than myosin V through its cooperative nature 
at zero load. These different mechanics are ideal for the respective myosin's physiological functions. 

PACS numbers: 87.16.Nn, 05.10.Gg, 82.39.-k, 82.37.-j 



Molecular motors play essential roles in various physi- 
ological functions, such as muscle contraction, molecular 
transport, cell motility and cell division through displace- 
ments powered by the chemical energy of ATP hydrolysis. 
The primary goal of molecular motor biophysical stud- 
ies is to understand how chemical energy is converted 
into mechanical movement, i.e. mechanochemical cou- 
pling. Single molecule techniques are a powerful tool 
to investigate chemical reactions and mechanical move- 
ments by molecular motors because they can directly ob- 
serve the elemental reaction process that is undetectable 
in ensemble average measurements P, Q ■ This is espe- 
cially true for motors like myosin V, which function as 
a single molecule or in cooperation with a small number 
of other myosin V molecules Q. By combining single 
molecule studies with in vitro kinetics studies, myosin V 
is thought to produce regular 36 nm steps coupled to one 
ATP hydrolysis ^4]. This means that the mechanochem- 
ical coupling of myosin V is a one-to-one relationship, 
i.e. "tight coupling", although no direct measurement 
has shown this relationship. In contrast, little is known 
about the mechanochemical coupling of ensemble func- 
tioning motors, such as myosin II in muscle. Because 
of its poor stepping ability, mechanical steps are eas- 
ily buried in thermal noise. This makes direct detec- 
tion of individual mechanical events difficult Q, leav- 
ing most to discuss only the average displacement Q. 
Furthermore, because in vivo myosin II work in ensem- 
ble, single molecule characteristics may not accurately re- 
flect the mechanochemical coupling of myosin II. In fact, 
multi-molecule systems show various characteristics un- 
predictable from single molecule studies [7, 8]. Therefore, 
the mechanochemical coupling of myosin II is still poorly 
understood. 

In order to examine mechanochemical coupling of both 
kinds of molecular motors, in this report, we investigated 
the characteristics of fluctuation in myosin V and myosin 
II motility and discuss their mechanochemical coupling 
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FIG. 1: Schematic of the experiments and typical traces of 
myosin motility, a, Single molecule motility assay of myosin 
V. The single molecule movement of myosin V along the actin 
filament was directly measured. Myosin molecules were genet- 
ically fused GFP (Green Fluorescent Protein) for observation 
(green star). Fluorescence from myosin in solution is not im- 
aged because of its fast diffusion whereas molecules moving 
along actin filaments are imaged as moving bright spots. The 
experiments were performed at 25 °C. b, In vitro actin gliding 
assay. The movement of an individual actin filament driven 
by a myosin II ensemble was measured. Actin filaments were 
attached to the fluorescent dye tetramethyl - rhodamine phal- 
loidin for observation (colored in orange). Myosin molecules 
were not labeled and not observed. The experiments were 
performed at 28 °C. c. Typical displacement time series ob- 
tained from the motility assays. ATP concentration, [ATP], 
was 100 fiM in both assays. 



consistently. Since these fluctuations contain the un- 
derlying molecular process, we can extract pertinent in- 
formation of mechanochemical coupling from fluctuation 
analysis, as Schnitzer et al. did to show one kinesin step 
couples to a single ATP turnover Q. 

In myosin V experiments, we directly observed sin- 
gle myosin molecule movement along actin filaments 
llOl . Fluorescence imaging was performed by TIRF mi- 
Actin filaments were adsorbed onto a 
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Freely floating fluorescently dyed myosins 
in solution bound to actin filaments commencing move- 
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ment. The movement was tracked by fitting the fluores- 
cent spots with a two-dimensional Gaussian distribution 
function. On the other hand, in myosin II experiments, 
we performed the in vitro actin ghding assay [l^ (Fig- 
[T|) and measured the shding velocity between the actin 
filament and myosin II molecules. Myosin II molecules 
were adsorbed onto the substrate. Introducing actin fila- 
ments with fluorescent dye and ATP resulted in a sliding 
movement between the myosin molecules and the actin 
filaments. The trajectories of the actin movement were 
obtained by tracking the center of the brightness of the 
actin fluorescence. For the convenience of the analy- 
sis, we chose the length of the actin filament to be less 
than 600 nm. In order to obtain velocity statistics, we 
converted two dimensional trajectories of both experi- 
ments into a displacement time series and calculated the 
Mean Square Displacement (MSD). We fitted the MSD 
to MSD(At) = ulAt'^ + alM + ^, where is the mean 
velocity, ct^ is the variance of velocity per unit time, ^ is 
the measurement error, and At is the sampling time. We 
performed both experiments at various ATP concentra- 
tions, calculated the MSDs for each [ATP], and identified 
the above parameter values. Here, we set Ai = 200 ms 
for myosin V and 32 ms for myosin II. For myosin V, 
MSDs were calculated from 100 traces, each of which 
contained more than 12 sampling points on average at 
[ATP] ranging from 2 to 1000 ^M; for myosin II, MSDs 
were calculated from 30 traces, each of which contained 
more than 100 sampling points on average at all [ATP] 
conditions. In very low [ATP] conditions (200 nM to 1 
/iM), we calculated myosin V MSDs from 40 traces, each 
of which contained more than 30 sampling points on av- 
erage at Ai = 2 s. 

The relationship between [ATP] and the mean velocity 
is shown in Fig. [2] . This relationship is known as the 
Michaelis - Menten (MM) mechanism, which represents 
the catalytic activity of the enzyme. We applied the MM 
mechanism to the actomyosin sliding movement : 

(Scheme 1) 

A-M + ATP ^ A-M ■ ATP A-M + Pr 

, where M, A, and Pr denote myosin, actin, and the ATP 
hydrolysis products ADP and Pi, respectively. In this 
scheme, an ATP binds reversibly to the actomyosin com- 
plex. The actomyosin catalyzes ATP hydrolysis to pro- 
duce a regular displacement, i.e. tight coupling. Note 
that in the case of myosin II, a single actin filament inter- 
acting with multiple myosin II molecules follows Scheme 
1, whereas in the case of myosin V, single myosin V 
molecule obeys MM kinetics. Although in the former 
case, the number of interacting myosin II molecules seems 
to affect the frequency of the chemical reaction, it was 
experimentally confirmed that the velocity of the single 
actin filament saturates and is independent of the num- 
ber of interacting myosin molecules except for myosin 
densities much low than those used here 12, 131 . In this 
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FIG. 2: The dependency of the mean velocity on the ATP 
concentration in motility assays. Solid lines show the fit of 
the Michaelis- Menten equation to the experimental data. See 
text for details. 



scheme, the ATP dependence of the sliding velocity is de- 
scribed by F = VraaAATP]/iKra + [ATP]) , whcrC Vma. 

is the maximal velocity at saturating [ATP] and Km is 
the [ATP] corresponding to the half maximal velocity 
of the sliding movement, called the Michaelis constant. 
Solid lines in Fig. [2] show the fitted MM equation with 
Vmax = 689 nm/s and Km = 24.6 /iM for myosin V 
and Vmax = 5960 nm/s and Km = 93.0 /iM for myosin 
II, respecti vely . These values are consistent with earlier 
studies [isLfl^. This shows that regarding mean velocity. 
Scheme 1 holds for both kinds of myosins. 

Then, is this simple MM mechanism also valid for the 
description of the velocity fiuctuation? In order to ex- 
amine this point, we further investigated the relationship 
between /i„ and a1 (Fig. [3]). In the case of myosin V, 
there is a linear relationship between /i^, and cr^ (Open 
circle) except for the velocity region corresponding to the 
[ATP] around the Michaelis constant. To test the valid- 
ity of Scheme 1, we performed stochastic simulations and 
compared these with the experimental results. Scheme 
1 requires three rate constants even though experimen- 
tally only two parameters, Vmax and Km, were obtained. 
Therefore we simplified Scheme 1 without losing any rel- 
evant statistics: 

(Scheme 2) 

A-M + ATP ^ A-M ■ ATP 

Since the left and right states in Scheme 1 are identical 
for the actomyosin complex, we can rewrite Scheme 1 as 
Scheme 2. Here we assumed a regular stepping motion 
is accompanied by a left directed reaction. Although the 
left directed arrow should include both ATP unbinding 
and ATP hydrolysis reaction, the former only reduces the 
stepping rate. This causes the reduction of the velocity at 
each [ATP], but does not affect the relationship between 
the mean and the variance of the velocity. Thus we can 
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FIG. 3: Relationship between and a'^. To show the depen- 
dence of the velocity variance per unit time on fj,v , we normah 
ized the following variables as follows : fiv = fJ-v / Vmax, cr^ 
= (Ti, / V^ax ■ Lines show the results of stochastic simulations 
of the MM mechanism. Parameters used in the simulations 
are as follows: fc_ = 18.5 and fc+ — 0.750 s~^/ii\f~^ 
for myosin V, fc_ = 596 and k+ = 6.40 s~^fj.M~^ for 
myosin II with 10 nm step, and fc_ = 30.0 and k+ = 
0.320 s"VA^"^ for myosin II with 200 nm. 

ignore the ATP unbinding reaction in this analysis. The 
rate constant of the left directed reaction, k- , was set to 
be fc_ =Vmax I step-size in order to obtain the exper- 
imental Vmax- The rate of the right directed reaction, 
A;+, was calculated by k^ = k- j K^- We calculated 
the number of stepping events for a unit time interval 
in the stochastic simulation with the Gillespie algorithm 
[ist . The step-size of myosin V was set at 36 nm, based 
on an earlier study 0]. The simulation results (red line 
in Fig. [3]) show good agreement with the experimental 
data. This confirms that myosin V is a tight coupling 
motor, consistent with earlier studies. 

The resulting relationship between /x„ and ti^ is ex- 
plained as follows : in very low [ATP] conditions, ATP 
binding is the single rate limiting step in the reaction 
sequence. Thus, the number of occurring chemical re- 
actions in a given time interval can be described by a 
Poisson distribution, in which the linear relationship be- 
tween and (7^ holds. On the other hand, when [ATP] is 
near the Michaelis constant, the frequency of ATP bind- 
ing and the hydrolysis reaction are comparable. Here, the 
fluctuation of each reaction cancels, and the fluctuation 
decreases. In higher [ATP] conditions, however, ATP hy- 
drolysis becomes a single rate limiting step and ct^ again 
becomes proportional to (We also confirmed the 

change of the number of the rate limiting steps through 
the randomness parameter analysis 0). 

In the case of myosin II, the relationship between 
and is qualitatively different from that of myosin V. 

is proportional to fiy in the lower velocity range, while 
ay is proportional to /Lt„ in the intermediate and higher 
velocity range. The crossover of the relationship occurs 



at the velocity corresponding to [ATP] = 2 /zM, which 
is similar to the Michaelis constant of ATPase [l3|. We 
also performed stochastic simulations for Scheme 2 for 
the myosin II data. Here, the movement of a single actin 
filament is assumed to be displaced 10 nm by individual 
myosin II molecules, as is expected from the tight cou- 
pling model [lij (Fig. [3] in green). The simulation results 
of scheme 2 showed a linear relationship between and 
l^v that could not reproduce the crossover or the larger 
fluctuations appearing in the experimental data. The 
change in the number of interacting myosin II molecules 
i.e. the number fluctuation cannot produce the crossover 
for the following reasons: the number fluctuation has a 
temporal correlation in the time scale of L/V, where L 
is the length of the actin filament and V is the velocity. 
This becomes smaller as the [ATP] increases. When con- 
sidering the velocity fluctuation per unit time, the effect 
of the number fluctuation also becomes smaller as the 
[ATP] increases because it is averaged out by the law of 
large numbers. Furthermore, to achieve the Vmax value 
with a 10 nm step, we have to set the stepping rate to 600 
s~^, which is large enough to average out the fluctuation 
of reactions. Thus, it is necessary to have a displacement 
larger than 10 nm per one ATP hydrolysis to reduce the 
number of stepping events per unit time and realize the 
larger fluctuations (for example, a 200 nm step simulation 
is shown in Fig. [3] in blue). This cannot be achieved in a 
tight coupling model because it is not physically plausible 
for the myosin molecule to generate a regular displace- 
ment larger than its own body («20 nm). Therefore, to 
produce a step-size more than 20 nm, multiple stepping 
events are required to occur in a single ATP hydrolysis. 

Employing a larger step-size can reproduce the large 
fluctuation as shown in Fig. [31 However, it still can- 
not reproduce the crossover relationship, because the 
stochastic property of Scheme 2 only holds the /x„ oc cr^ 
relationship, as mentioned above. Therefore, we incorpo- 
rated a variable step-size depending on the [ATP] , taking 
into account the change in the relationship between /i„ 
and al around the Michaelis constant of the ATPase, 
2 /iM. Here we assume that the ATPase kinetics fol- 
lows the Michaelis constant of 2 /iM. This assumption 
leads to the relationship Vmot — VATPase x d, where 
d is the step-size per one ATP hydrolysis. Following 
from Scheme 2, Vmot = yr,°i[ATP]/(if ™°* + [ATP]) 
and Vatp = V^IP[ATP]/{K^^p +[ATP]), then we ob- 
tain d = dmaxiK^^'' + [ATP])/{KZ°' + [ATP]), where 
dmax — V^°llV^^ . In this expression, d is variable 
depending on the [ATP]. In fig. IH we show stochas- 
tic simulations using scheme 2 with a step-size function. 
Rate constants were set as in the case of Fig. [Sj except 
for fc+ = A;_ / K^"^^ . We employed dmax to be 100, 
200 and 300 nm to compare with the experimental data, 
making dmax five to fifteen times larger than the size of a 
myosin molecule. The large dmax is expected to be due to 
stochastic multiple steps by myosin molecule(s), so it is 
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actually more plausible that dmax distributes with a large 
mean value (100 '~ 300 nm) rather than being constant. 
FigU] shows a simulation with an exponential distribu- 
tion and a mean of dmax- A large but constant d^ax 
(200 nm) can also explain the dependency of cr^ on /i^, 
but not as well as the exponential distribution assumed 
here. Overall, the result shows that the model could suc- 
cessfully reproduce the experimental data with the most 
likely maximum step-size, dmax, being 200 nm. Thus, the 
characteristic the step size dependence on the [ATP] is 
essential to reproduce the crossover relation, confirming 
the "loose coupling" nature of myosin II [17]. Note that 
it is not essential assumption that the step is exponen- 
tially distributed, because the simulation result in which 
the step size is exponentially distributed but not depen- 
dent on the [ATP] condition does not show the switching 
behavior (Fig. IH in purple). Here, we should again stress 
that our analysis does not need to estimate the number 
of interacting myosin molecules to obtain the displace- 
ment of an actin filament per ATP hydrolysis cycle. Al- 
though earlier studies reported the step-size, they needed 
to estimate the number of interacting molecules and this 
discrepancy lead to the one-order difference in calculated 



step-size [12, ll3j . Estimating the number of interacting 
molecules is extremely difficult and therefore unreliable. 
Our method can overcome this difficulty, and clarify the 
manner of the mechanochemical coupling. This confirms 
that our fluctuation analysis is a promising tool to inves- 
tigate the mechanochemical coupling in molecular mo- 
tors. Also, we should note that dmax ~ 200 nm is not 
unrealistic because Yanagida et al. and Higuchi et al. 
[isl . [iqI . [20| have already suggested the likelihood of a 
value much greater than that predicted from the tight 
coupling model fl^ . 

The resulting characteristics for myosin-V and myosin- 
II are reasonable considering their physiological func- 
tions. Myosin-V transports cargo over a long distance 
by a single or a small number of motors. The tight 
coupling mechanism guarantees a consistent step-by-step 
movement over this long distance. On the other hand, 
myosin-II works in an ensemble like muscle to produce 
adaptive motions depending on the condition. The loose 
coupling mechanism is advantageous for such adaptation. 
For example, this mechanism can generate a large force 
with a short step at a high load but also a large step at a 
low load to conserve energy. Because a single rnyosin-II 
motor cannot generate steps larger than 30 nm [5| , steps 
as large as 200nm should be produced not by a single 
myosin-II motor but by the cooperation of many [^l'] . It 
is possible that the active steps by some myosin molecules 
cause passive steps in others through the backbone con- 
necting motors. To elucidate how the cooperative dy- 
namics arise from the interaction of multiple motors and 
how this adapts to the biological environment is the focus 
our future work. 
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FIG. 4: Comparison with our proposed model to the experi- 
mental data of myosin II. Symbols are the same as the myosin 
II data in Fig. O Solid lines show the simulated relationship 
of our model. See text for details. Inset, The dependence 
of the mean velocity on the [ATP] . Symbols are the same as 
myosin II data in Fig. [2] Red solid line is the stochastic sim- 
ulation. The parameters are dmax = 200 nm, — 2fiM, 
= 93^M, respectively. 
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